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Abstract 

We investigate hard radiation emission in small-angle transplanckian scattering. We show 
how to reduce this problem to a quantum field theory computation in a classical background 
(gravitational shock wave). In momentum space, the formalism is similar to the flat-space 
light cone perturbation theory, with shock wave crossing vertices added. In the impact 
parameter representation, the radiating particle splits into a multi-particle virtual state, 
whose wavefunction is then multiplied by individual eikonal factors. 

As a phenomenological application, we study QCD radiation in transplanckian collisions 
of TeV-scale gravity models. We derive the distribution of initial state radiation gluons, 
and find a suppression at large transverse momenta with respect to the standard QCD 
result. This is due to rescattering events, in which the quark and the emitted gluon scatter 
coherently. Interestingly, the suppression factor depends on the number of extra dimensions 
and provides a new experimental handle to measure this number. We evaluate the leading- 
log corrections to partonic cross-sections due to the initial state radiation, and prove that 
they can be absorbed into the hadronic PDF. The factorization scale should then be chosen 
in agreement with an earlier proposal of Emparan, Masip, and Rattazzi. 

In the future, our methods can be applied to the gravitational radiation in transplanckian 
scattering, where they can go beyond the existing approaches limited to the soft radiation 
case. 
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1 Introduction 



Scattering at center-of-mass (CM) energies exceeding the quantum gravity scale (transplanckian 
scattering, or T- scattering, for short) is an exotic process of significant theoretical interest. In 
particular, it provides a laboratory to study the black hole information loss paradox. Microscopic 
black hole formation and its subsequent evaporation is expected for impact parameters b of the 
order of the Schwarzschild radius Rs of a black hole of mass y/s [I] , [2] , [3] , |1] . The detailed 
description of how this happens depends on the unknown underlying theory of quantum gravity 
and is at present out of reach. On the other hand, large impact parameters b ^ Rs correspond to 
elastic small-angle scattering, whose amplitude can be predicted on the basis of General Relativity 
alone. It is given by eikonalized single-graviton exchange [I] , [5] , [6] , [7j . Computing the corrections 
in b/Rs to the elastic scattering, one hopes to learn about the strong inelastic dynamics at 6 ~ Rs 

IH],i],[in!,III],II3. 

T-scattering is also interesting phenomenologically. If large extra dimension scenarios of TeV- 
scale gravity [13] are realized in Nature, this process could be observed at the LHC and other 
future colliders [ll],[l5], as well as in collisions of high-energy cosmic neutrinos with atmospheric 
nucleons [TB] , [Uj . In these scenarios the total T-scattering cross section is finite, grows with energy, 
and is dominated by calculable small-angle scattering between partonic constituents [IT], [15], see 
Fig. [Tj The subleading black hole production cross section at present can only be estimated from 
geometrical arguments. 





Figure 1: The up T-scattering in cosmic ray physics (left) and the pp T-scattering at 
the LHC (right). In both cases the dominant process is small-angle elastic scattering 
between partons, giving jet(s)-\- anything in the final state. The zigzag line denotes the 
eikonal 2^2 scattering amplitude. 



In spite of the small scattering angle, the typical momentum transfer in these scattering events 
is well above the QCD scale, and the typical impact parameter is much smaller than the proton 
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size, which sets the typical distance between two uncorrelated partons inside the proton. It is 
unhkely that a multiple parton interaction, Fig. |2| will occur in the same T-scattering. Thus it 
is clear that the partonic picture should be applicable at leading order in the QCD coupling. In 
other words, we can compute the total cross section via a convolution of the partonic cross section 
and the parton distribution functions (PDF) f{x,fip). 




V 



Figure 2: Multiple partonic interactions, in which more than one pair of partons ex- 
change gravitons, are suppressed. This is because the eikonal phase vanishes very 
quickly with the transverse separation, and it is unlikely to find a second pair of partons 
for which this phase is large. 

Several interesting questions arise when one tries to think what happens beyond the leading 
order. For example, it's not known how to treat events of the type shown in Fig. [3} where one 
of the colliding partons (say a quark) radiates a gluon just before the collision. Now the quark- 
gluon separation is not necessarily large, and the pair may scatter coherently. What is the correct 
description of such rescattering processes, and what is the resulting effect on the total T-scattering 
cross section? 

A related question is: which factorization scale fip should we choose in the computation 
of the total cross section? As is well known, QCD-initiated processes have significant higher- 
order logarithmic corrections, associated with the coUinear QCD radiation off initial partons. By 
choosing ij,p appropriately, these corrections can and should be reabsorbed into the PDF. As we 
will see below, the familiar choice fi'jp ~ —t is likely not the right one for the T-scattering. If so, 
we would like to see this explicitly. 

The purpose of this paper is to answer the above-mentioned questions. We will be focusing on 
the QCD radiation since it is the dominant phenomenological effect due to the relative largeness 
of as- However, our methods are equally applicable to the radiation of photons or any other spin 
1 gauge bosons. We also hope that these methods may later prove useful in the more complicated 
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Figure 3: Initial state radiation processes, when a quark emits a gluon just before par- 
ticipating in T-scattering, are expected to play a role. They give logarithmic corrections 
to the cross section, which determine the optimal factorization scale fip of the process. 

problem of gravitational radiation emitted in T-scattering, and in particular to provide an alter- 
native to the existing computations which are limited to the case when the emitted radiation is 
soft ^M- 

The paper is organized as follows. In Section [2] we review the eikonalization of small- angle 2 
2 partonic scattering amplitude, following [17], [15]. To keep close contact with phenomenology, 
we work in the context of large extra dimension scenarios with the quantum gravity scale around 
a TeV. A key feature of the extra-dimensional situation with n compactified dimensions is the 
appearance of a new length scale be, which sets the range of a typical T-scattering interaction. In 
D- dimensional Planck units, D = 4 -|- n, we have 



6c ~ s 



i/ji 



so that be ^ Rs in the deep transplanckian regime a/s ^ 1. We also present an alternative 
computation of the 2^2 amplitude, based on generalizing to D dimensions the early idea of 't 
Hooft [T], who considered the small-angle T-scattering by solving the Klein-Gordon equation for 
one particle propagating in the classical gravitational field of the other particle. 

In Section |3] we start discussing small-angle T-scattering with hadronic initial states, as in 
Fig. [T| in which QCD effects are expected to play a role. According to the existing proposal of 
Emparan, Masip and Rattazzi [T7], the total cross section for these processes must be computed 
with the following prescription for the PDF scale fip '■ 



M(i) = \ ^,,1,^ ,_i (1-1) 

[ [b^ ^)"+^ if q> b^, 
where q = a/— For sufficiently high momentum transfers this deviates from the familiar pre- 
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scription //^ ~ g. An intuitive justification for tliis scale in terms of the typical impact parameter 
was given in [TTj, but we would like to check it via a direct computation. 

The first step is to be able to evaluate the amplitude for small-angle T-scattering accompanied 
by coUinear QCD radiation. In the resummation approach [I7],[l5], this computation seems 
prohibitively difficult even for one-gluon emission. Indeed, the eikonal amplitude for quark-quark 
T-scattering is a sum of an infinite number of crossed ladder gravition exchanges. The outgoing 
gluon may be attached anywhere on the quark lines, both external and internal. Moreover, in the 
q q + g splitting, the emitted near-collinear gluon is not necessarily soft, and thus may also 
exchange gravitons. The number of diagrams to resum skyrockets. 

't Hooft's approach is a much better starting point. As we point out, it can be easily 'upgraded' 
to the case when radiation is present, provided that only one of the two colliding particles radiates. 
This covers completely lepton-quark scattering and is an important special case for quark-quark 
scattering. The idea is very simple. In the 2 — 2 scattering, 't Hooft treated one particle 
classically, the other one quantum mechanically. The only new twist is to allow the quantum 
particle to radiate. In other words, we should treat the non-radiating parton classically, while the 
radiating parton and the gluonic radiation field with which it interacts quantum mechanically. 

This trick reduces the problem to a quantum field theory computation in the classical gravita- 
tional background produced by a relativistic point particle, the Aichelburg-Sexl (AS) shock wave 
jl9j . In Section |4] we develop the necessary formalism. We first consider the simplest perturbative 
quantum field theory in the AS background: a scalar field with cubic self-interactions. We in- 
troduce a diagram technique for computing arbitrary transition amplitudes in this theory, which 
turns out to be closely related to the standard rules of light-cone perturbation theory in flat space. 
We then explain the changes necessary for the gluon field and for the scalar-gluon interactions, 
and compute the one-gluon emission amplitude as an example. 

Notice that while fermionic matter fields can be considered analogously, we do not include 
them in this work in order to keep technical details to a minimum. Thus we stick to a toy model 
in which the partonic constituents of colliding hadron(s) are scalars. 

Armed with the knowledge of gluon emission amplitudes, in Section |5] we attack the question 
of QCD corrections to T-scattering. For definiteness and simplicity, we consider the gravitational 
analogue of the DIS: a transplanckian electron-proton collision. The observable is the total cross 
section as a function of the Bjorken x and the transverse momentum transfer q, in the small- 
angle region |q| <^ y/s. At leading order (LO) in the QCD coupling a^, the partons scatter 
elastically on the electron (no gluon emission). At next-to-leading order (NLO), we demonstrate 
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the appearance of logarithmic corrections whose scale is precisely the /i_p(i?) from Eq. (1.1). We 
find that the cross section factorizes, in the sense that these logarithms appear multiplied by the 
DGLAP splitting functions, and can be reabsorbed into the PDFs. Finally, we are able to show 
that this factorization holds to all orders in a^, in the leading-logarithm approximation (LLA). 

Our computation gives an explicit check for the validity of the partonic picture for the T- 
scattering. Moreover, it gives an interesting and unexpected explanation for why the PDF scale 
deviates from the usual /ii? ~ g. It turns out that rescatterings like in Fig. |3] suppress the initial 
state QCD radiation at transverse momenta q > b^^. As a result the transverse momentum 
distribution of emitted gluons has the form: 

where {dN / d'^q)^ ~ 1/q^ is the standard distribution without rescattering, and /(g) is a function 
interpolating between 1 for g ^ h^^ and l/{n + 1) for g ^ h^^. Logarithmic corrections to the 
cross section are obtained, as usual, by integrating Eq. (1.2) over the gluon phase space, and the 



scale of these logarithms is a geometric mean of g and as in Eq. (1.1). 

Notice that one could imagine other distributions giving rise to log/ii?(g), for example the 
standard 1/q^ with a sharp cutoff at h^^ . In this sense Eq. (1.2) contains more information 
than the identification of the correct factorization scale. The predicted suppression of the initial 
state radiation is n-dependent and could in principle be used to determine the number of extra 
dimensions. 

A crucial insight into the physics of radiative processes is obtained by going into the im- 
pact parameter representation. In this picture, we find that the scattering is described via a 
multi-particle wavefunction of the virtual state (parton+radiated quanta), which is multipled by 
individual eikonal factors when crossing the shock wave. This interpretation suggests a possible 
generalization of our formalism to the case when both colliding partons radiate, which we discuss 
in Section [6l 

In conclusion, this work shows that factorization holds for QCD effects in T-scattering, and that 
the factorization scale has a nontrivial dependence on g^, in agreement with the earlier proposal 
of Ref. llTj. The novelty is that we arrive at these results by a concrete computation, and that 
we derive the modified distribution of the initial state radiation due to rescattering effects. The 
new distribution should be now incorporated into a 'transplanckian parton shower algorithm', to 
be used in Monte-Carlo simulations of T-scattering. We will come back to this issue in a future 
publication. 
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2 Review of the eikonal approach to T-scattering 



In this section we review the basics of small angle T-scattering in the eikonal approximation. We 
will work within the large flat extra dimensions scenario of TeV-scale gravity [13], see [20] for the 
current experimental constraints. 

Consider then two transplanckian massless Standard Model (SM) particles, thus confined to 
the SM 3-brane, which scatter due to the Z)-dimensional gravitational field, D = A + n, n being 
the number of large extra dimensions, n > 2 for phenomenological reasons. For now we ignore 
all interactions except for gravity. In particular, we suppose that the colliding particles are not 
charged, and thus cannot emit photons or gluons. We are interested in the scattering amplitude 
for small momentum transfer —t/s <^ 1. In this regime gravitational radiation is also suppressed 
(see [15]), and we have elastic 2^2 scattering. 



2.1 Resummation 

The most direct way to compute the amplitude is by resumming the crossed ladder graviton 
exchange diagrams [I7],[l5], see Fig. |4j For small momentum transfer, exchanged gravitons are 
soft, and well-known simplifications occur in the vertices and the intermediate state propagators, 
allowing the resummation. The first term in the series, the one-graviton exchange, is given by 

s2 r 



^Born(q) = -^y (2.1) 

where q is the momentum transfer, which lies mostly in the direction transverse to the beam: 
t ~ — q^- The D-dimensional Planck scale M^) ~ 1 TeV is normalized as in [I5],[20]. The divergent 
integral over the extra dimensional momentum / needs to be treated properly; see below. The 
second term in the series, the sum of two one-loop diagrams, turns out to be equal to a convolution 
of two Born amplitudes: 

^l-loop(q) = ^ y ^^;^^Born(k)v4Born(q " k) , 

and this pattern continues to higher orders. As a result the series can be summed by going to the 
impact parameter representation. The amplitude acquires the eikonal form: 

Aik(q) = ^Born + ^1-ioop + • • • = -2is j d% e-'^"^ {e'^ - 1) , (2.2) 

with the eikonal phase x given by the Fourier transform of the Born amplitude in the transverse 
plane: 
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Figure 4: The 2 — > 2 small angle T-scattering amplitude is given by a sum of crossed- 
ladder graviton exchanges. 



To evaluate the eikonal phase, we need to regulate the divergent Born amplitude (2.1 ). In 



dimensional regularization was used, and it was argued that since the subtracted divergent terms 
are local, they do not affect the small angle scattering amplitud^ The eikonal phase was found 
to be: 



X{h) 



1 



D 



(47r)t-T(n/2) 



l/n 



Ml 



l/n 



(2.3) 



The corresponding amplitude is then given by: 



Fniy) = -I 



ATishlFn{hM\) , 
dx xJo{xy)[e^^ 



(2.4) 



The functions Fn{y) are plotted in Fig. [5] (see also Fig. 2 of [15J). Their most salient features are 
as follows. At moderate y < 1, we have Fniy) = C>{l),^the integral (2.4) receiving contributions 
from x ~ 1. On the other hand, for y ^ 1 the integral has a saddle point at a;* = (n/?/")"+i <^ 1, 
and the amplitude decays: 

1 

, ^ , , , 72 " + 1 _n + 2 

\Fniy)\^^==y (i/>i)- 

The appearance of the scale be is a peculiar feature of T-scattering for n > 0. Since the 
amplitude is the largest in the region y < 1, a. typical scattering will have |q| < b^^. Yet a 
classical particle trajectory for these q is undefined, all impact parameters b be contributing 
to the scattering. On the other hand, for 6c|q| ^ 1 (y ^ 1) the scattering is dominated by 
a characteristic impact parameter b^, = bcX^,, corresponding to the above saddle point. In this 

"'^This was later confirmed in |21j by using in (2.11 a physical regulator exp(— fc^w^), with w interpreted as an 
effective width of the SM brane, w ^TeV^"'^. It was found that the resulting eikonal phase coincides with (2.3 1 for 



b > w, while for b < w it varies slowly (logarithmically). The eikonal amplitude then indeed agrees with (2.4) in 

the small scattering angle region. The same conclusion was also reached in [22| using a sharp cutoff. 
^For n — 2 there is mild logarithmic growth. 
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Figure 5: The functions Fn{y) for n = 2,3,4,5; absolute value (solid), real part 
(dashed), imaginary part (dot-dashed). Notice that QmFn{0) < oo for all n > 2, 
implying finite total cross section. On the other hand, 3?ei^„(0) < oo for all n > 3. 
See JT^. 

case the particle trajectory is well defined and the T-scattering is truly semiclassical, with many 
gravitons being exchanged. 

2.2 't Hooft's method 

An alternative computation of the small angle T-scattering amplitude can be given using a method 
due to 't Hooft [1] , originally formulated in four dimensions. In this approach, particle B scatters 
on the classical gravitational field created by particle A. In other words, particle A is treated as 
a classical point particle, while particle B is treated quantum-mechanically. 

Consider then the gravitational field of a relativistic classical point particle A of energy Ea 
propagating in the positive z direction. This field is the D-dimensional generalization of the AS 
|19j shock wave: 

ds'^ = -dx~^dx~ + ^{x±)S{x~){dx~y + dx]_ . (2.5) 
Here x"^ = t ± z, while x± denotes D — 2 transverse directions. Einstein's equations with the 
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light like source 

reduce to one linear equation for the shock wave profile $: 

- 91$ = 167rGDEA5^^'^\x±) . 



(2.6) 



The solution of this equation coincides with the eikonal phase (2.3) per unit of particle B energy: 



$(xx) = Ej^\{x±) . 



(2.7) 



The right-moving particle B is confined to the SM 3-brane, and its wavefunction solves the 



Klein-Gordon equation in the metric induced on the brane by the shock wave (2.5). At a: < 
the wavefunction is a standard plane wave 

0(x) = exp{ipB-x) = exp{—iEBX~^) . 



The metric (2.5) has a strong discontinuity at x =0. To solve the Klein-Gordon equation 



across the discontinuity, it is convenient to make a coordinate transformation [3].[23] 



X 



X 



X 



x+ + e{x-)^ix,) + x-0(x-)^^^^^ , 



(2.8) 



Xi 



In the X coordinates the metric is continuous across x~ = 0. When crossing the shock wave, the 
wavefunction remains continuous in these coordinates. This means that for small positive x~ we 
have: 

0(x) = exp{-iEBX+) = exp[-iEB{x^ - $(x)] . (2.9) 

The x-dependent shift of the x^ coordinate has a well-known classical origin: it is related to the 
time delay experienced by geodesies crossing the AS shock wave, see Fig. pf\ 



We now see from (2.9) that the wavefunction immediately before and after the collision is 



related by a pure phase factor exp(ii?B$(x)), which via (2.7) is identical with the eikonal phase 



^We would like to stress the auxiliary character of the x coordinates, in which the metric is not manifestly 
flat for a;^ > 0, nor even manifestly asymptotically flat, which makes these coordinates unsuitable for deflning 
asymptotic outgoing states. The asymptotic states should be described in the x coordinates, that's why in the last 



equation in (2.9) we reverted to them. 
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Figure 6: The gravitational field of particle A (the AS shock wave) is concentrated on 
the null plane x~ = 0. Geodesies crossing this field experience an ^-dependent shift of 



the x"*" coordinate. This is the same shift as in Eq. (2.9) 



factor in (2.3). An alternative derivation, by directly solving the Klein-Gordon equation, is given 
in Appendix |Aj 

Thus, 't Hooft's method is equivalent to the resummation. This is not surprising, because the 
external field approximation in quantum field theory resums precisely crossed ladder diagrams 
|24j . The AS shock wave is a solution to both linearized gravity and the full nonlinear Einstein's 
equations. In retrospect, this explains why the diagrams in which gravitons emitted by particles 
A and B interact did not have to be taken into account in the resummation method. See [7] for 
a detailed discussion and comparison of the two methods m D = A. 

Still, an attentive reader will notice two small differences between the two results. First, 



Eq. (2.3) contains —1 under the integral sign, while 't Hooft's method gives a pure phase. This 



is the usual difference between the S- and T-matrices, 5 = 1 + iT. Second, the amplitude (2.3) is 
relativistically normalized, while in the new derivation normalization needs yet to be determined. 
Modulo the normalization issue (which will be resolved in Section 4. 1| below), the power of 't 



Hooft's method relative to the resummation is quite evident. The eikonal phase is given a simple 
physical interpretation — it is related to the time delay experienced by geodesies upon crossing the 
shock wave. The exponential factor e^^ emerges as a whole rather than by summing infinitely 
many individually large terms. 



3 T-scattering with hadrons: intuition and questions 

If TeV-scale gravity is the way of Nature, then transplanckian collisions may be within the en- 
ergy reach of the LHC Moreover, transplanckian collisions may be constantly happening in the 
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atmosphere, between the atmospheric nucleons and high-energy cosmic rays {^/s ~ 10^ GeV for 
Ecr ~ 10^^ GeV of the order of the GZK cutoff). In case of cosmic ray neutrinos this signal could 
actually be observable. 

Since protons are not elementary particles, the theory of small angle T-scattering from Section 
|2] should be applied instead to 2 ^ 2 collisions between the partonic constituents. Notice that 
since we are dealing with CM energies well over a TeV, the typical momentum transfers will be 
hard compared to the QCD scale, even though the scattering angle has to be small for the eikonal 
approximation to be valid. Thus the collision resolves the internal structure of the proton(s), and 
the partonic picture is applicable [17]. 

Viewed another way, when two protons collide, there is a phase factor ~ $(x — y) for each 
pair of partons moving in the opposite directions, see Fig. |2} This factor tends to zero rapidly at 
transverse separations |x — y| ^ be, where be ~ (100 GeV)~^ for T-scattering at the LHC energies. 
Since partons are distributed in the disk of radius (GeV)~^ ^ b^, it is unlikely that more than 
one pair will undergo a hard collision. 

We would like to briefly mention which observables one usually computes in phenomenological 
studies. In up collisions one is mostly interested in the total interaction cross section as a function 
of the energy transfer to the proton [H] , [22] , • We will discuss a similar observable in Section [5] 
below. On the other hand, in the pp collisions at the LHC one studies two jet final states of high 
invariant mass, produced at a small angle to the beam [I5],[26], see Fig.[T| These jets originate from 
all possible parton pairs {qq, qg, gg) with the same partonic cross section, the eikonal amplitude 
being independent of the particle spin. For M^, not much above a TeV, the dijet T-scattering 
signal turns out to be visible over the QCD background. 

So far it may look that from the point of view of QCD, the T-scattering is just like any other 
hard process. Let us however discuss which parton distribution factorization scale fip one should 
use when evaluating the T-scattering cross sections — a necessary prerequisite for any practical 
computation. 

For the usual hard processes, we are accustomed to the choice fip ~ |q|, but for the T- 



scattering this turns out to be more subtle. As we discussed in Section |2.1[ T-scattering becomes 
semi classical in the region of large momentum transfers |q| ^6^^. In this regime, the transverse 
distance characterizing the process is the typical impact parameter b^ ~ ^c/(^c|q|) which is 
parametrically larger than |q|~^. It is for this reason that Ref. [l7j advocated a hybrid prescription: 
one should use /i^ ~ |q| for |q| < b^^ and switch to jj,F ~ b^^ for |q| ^ 6^^, see Eq. ( l.l[ )j^ 



*A limiting case of this prescription in the context of black hole production was advocated earlier in the first 
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For T-scattering at the LHC energies, the factorization scale Hf{<i) will be hard with respect 
to the QCD scale as long as the momentum transfer q is hard. This gives a self-consistency check 
on the proposed picture]^ 

The above is a summary of the current understanding of QCD effects in T-scattering. Clearly, 
it is based mostly on intuition. We would like to develop a systematic theory of these phenomena. 
In particular, such a theory should allow to check the factorization scale proposal by a concrete 
computation. We have to evaluate the leading log corrections to the T-scattering cross section 
due to the initial state radiation emission, and to show that they can be absorbed into a shift of 
the PDF factorization scale. Since fip is conjectured to have a nontrivial dependence on q, some 
nontrivial physics is likely to come out. 

Two equivalent methods were given in Section [2] to describe T-scattering without radiation. 
Which one shall we try to generalize to the case when radiation is present? 

For the resummation method, generalization does not seem to be easy, not even for the one- 
gluon emission. Think about infinitely many crossed-ladder diagrams, infinitely many places to 
attach the gluon line, and the necessity to take into account the gravitational exchanges of the 
emitted gluon! 

For 't Hooft's method, on the other hand, the situation looks hopeful: if only particle B radi- 
ates, it is quite clear how to include its radiation. Namely, we should keep working in the classical 
gravitational background created by particle A, but switch from relativistic quantum mechanics 
(wavef unctions, the Klein-Gordon equation) to quantum field theory (Green's functions and in- 
teraction vertices). We will follow this path and will see that it allows relatively straightforward 
computations of the gluon emission amplitudes. 

Physically, the assumption that particle A does not QCD-radiate is realized if A is a lepton. 
If both A and B are strongly interacting, one could first compute the radiation off B (taking A 
classical), then off A (taking B classical). Such an approximation of independent emission is valid 
for the dominant, collinear, radiation in the usual perturbative processes. For the T-scattering, 
we will be able to partially justify it below. But first we have to understand well the case of 
non-radiating A. 

Ref. [H]. 

^Notice however that (g) is a decreasing function of the CM energy, since the typical impact parameter grows 
with ^/s. Even for hard momentum transfers, for sufficiently high ^/s the factorization scale will come down to a 
GeV, signalling a breakdown of the partonic picture. At even higher CM energies (which are well beyond the range 
of LHC or even cosmic ray collisions), the proton should interact gravitationally as a point particle. 
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4 Quantum field theory in the shock wave background 



4.1 Scalar field 

To compute the QCD radiation accompanying a transplanckian collision, we will replace particle 
A with the classical background it generates, but will keep particle B and the gluons as quantum 
fields. Thus we will be doing perturbative QFT computations in the shock wave background. We 
start with the simplest interacting QFT, the massless 0^ theory: 

C = y^{lg^''d,<pd.<P-^<p'). (4.1) 

We will describe how to compute transition amplitudes in this theory, and how these are related 
to the amplitudes in the full theory (i.e. before particle A was replaced by a classical gravitational 
field). 



The g^^ in (4.1) is the 4-dimensional metric obtained by restricting the D-dimensional AS 



shock wave (2.5) to the SM brane on which both particles and the radiation propagate. We will 



continue using the coordinates as in (2.5), only restricting the number of x± components from 
D — 2 to 2. Two features make this theory much simpler than it would be for generic curved 
backgrounds treated in [29j: 

1. the metric is invariant under shifts. The conjugate momentum p~ is conserved. This 
leads in particular to the absence of spontaneous particle creation. 

2. the spacetime is flat except on the = plane. The Feynman rules are simplifled by using 
the flat-space coordinates. 

We start by canonically quantizing the quadratic part of the lagrangian. The scalar fleld modes 
are found by solving the equations of motion (EOM) in the shock wave background with the plane 
wave conditions in the asymptotic past x~ < 0|^ 

^i"^ (x) = ^(-x-)e*[^'l-^ + e{x-) [ I{p-, q) e*[P+'i]-^ , (4.2) 



{2n) 



The compact "vector in square brackets" notation denotes an on-shell 4-vector whose + component 
is computed in terms of the known — and ±, i.e. [p + q] = ((p + q)^/p^,p^, p + q), etc. 



^By boldface letters p,x we denote the 2-dimensional, transverse to the beam, part of 4-diniensional Lorentz 
vectors. The Minkowski space signature is — h + +. 
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The function I{p , q) is identical to the eikonal amphtude (2.3), up to the normahzation and 



the absence of —1 under the integral sign (which means that it contains an extra 5-function piece). 



The modes (4.2) solve the Klein-Gordon equation both for x < and for x > 0. Across the 



shock wave, they satisfy the matching condition of Section |2.2 

<^(a:~ = +e, x"*", x) = </)(x~ = —e, — $(x), x) 
We proceed to quantize the field by expanding in oscillators: 

J 



dp d'^p 



Pi Pi' P2P2' 



Such normalization of the creation/annihilation operators is standard for quantizing on the light 
cone; it differs from the usual one by a simple rescaling. 

Equivalently, we can quantize using the outgoing modes, which reduce to plane waves for 
X' > 0: 



(x) = e{x~)e'^P^-'' + 9{-x- 



The in and out modes are related by a unitary Bogoliubov transformation, which acts only on 
the transverse momentum p but not on p~ . Thus there is no spontaneous particle creation in this 
background; the vacuum is unambiguously defined. 

Let us now build a perturbation theory for transition amplitudes. The logic is simplest in the 
position space. Even though the metric is singular at x~ = 0, it is easy to see thaty^ = 1: the 
metric determinant drops out of the interaction lagrangian. Thus the Feynman diagrams will be 
given by flat space integrals, with no singular contribution from the shock wave. For instance, the 
t-channel diagram contributing to the pi,p2 PsjPa, transition amplitude will be given by: 

-^^)' [d'xd'y[<f>'^-\x)Y[(j>T\y)rG{x,y)<pf{x)<jy^{y). (4.4) 



The 0°"* and 0"^ enter as the in and out state wavefunctions. The propagator G{x,y) must be 
a;~-ordered: 

Gix, y) = e{x- - y-) (O|0(x)</.(y)|O) + {x ^ y) . 
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In (4.4) we have to integrate in all possible x~ orderings of x and y with respect to each other 
and to the shock wave sitting at x~ = 0. The propagator will take different forms depending on 
the ordering. For x~ and y~ on the same side of the shock, we get the flat space result: 

'p->0 



Gflat(x,i/) = e{x- >y->0) I 2p-{2nY ^^'^^~'^ + v) ■ 



On the other hand, across the shock wave we have 

G,,U^,y) = 9ix- >0>y-) >o / (2^ ^(P^' + ^) • 

The momentum-space Feynman rule can now be found by straightforward Fourier transforma- 
tion; they are as follows. The in and out states are specified by the p~ and p of all incoming and 
outgoing particles. The transition amplitude i —>■ f in the external gravitational field of particle 
A is then given by: 

o„t {m,, = 2i2n)5ip]-pr)Mit^f). 

The f) is a function of the external momenta computed as a series in A according to the 

following rules. To obtain the O(A^) term: 

• Draw the Feynman diagrams with (j)^ vertices, considering all possible a;~-orderings of 
these vertices with respect to each other and to the shock wave at x~ = 0. 

• Consider all shock wave crossings as additional vertices, with entering transverse momenta 

representing momentum exchange with the shock wave. 

• Assign p~ ,p internal lines momenta by using their conservation in all vertices (0^ and shock 
wave crossings). Momentum flow is in the direction of increasing x~ . The internal p'^ 
momenta are not conserved but are assigned by using the on-shell condition p+ = p^/p~. 

• For each (j)^ vertex multiply by —iX. 

• For each shock wave crossing vertex multiply by p~I{p^,cia), where is conserved in the 
crossing. 

• For each internal line (i.e. a line connecting two vertices, 0^ or shock wave crossing) carrying 
momentum p~ , multiply by 0{p~)/p^ . 
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• The 0^ vertices and the shock wave at x~ = divide the x~ axis into two unbounded and 

bounded intervals. For each bounded interval, we define an intermediate state, consisting 
of all the particles whose internal lines traverse this interval. For each intermediate state at 
negative x~ , the amplitude is multiphed by 

i 

For each intermediate state at positive x~ , it is multiplied by 

i 

-EintermP+ + ^^ ' 

The sums are over all particles in the initial {x~ = — oo), intermediate, and final {x~ = +00) 
state. 

• Integrate over the momenta exchanged with the shock wave: 

/(2,)^i<^'(^q„ + p.-p,)ng;^. 

• For loop diagrams, integrate over all undetermined momenta (/c~,k): 

J W^' 

The reader will notice a striking similarity to the usual light-cone perturbation theory (PT) 
rules [30]. Notice in particular the light-cone energy denominators, and the 0{p~) factors, which 
eliminate some of the diagrams present in the time-ordered 'old' perturbation theory. New fea- 
tures in our case are the shock wave crossing vertices, and that there are two types of energy 
denominators, depending on the ordering with respect to the shock wave. We thus have 'light- 
cone PT in presence of an instantaneous interaction'. Many years ago, Bjorken, Kogut and Soper 
j31j have developed light-cone PT in external electromagnetic field, and argued that at sufficiently 
high energies interaction with the external field can be represented as an instantaneous eikonal 
scattering]^ In our case, the eikonal factor has gravitational origin, but the formalism is the same. 
The formalism of pT] has found application in the dipole scattering approach to the DIS at small 
x: an almost-real photon splits into two quarks which then undergo eikonal scattering in the gluon 
field of the proton [32]. The difference is that the gluon field of the proton is not really known, 
while in our case the eikonal phase can be computed exactly. 
^We thank Zoltan Kunszt for bringing this work to our attention. 
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We will now demonstrate the rules by computing a couple of amplitudes. The elastic one- 
particle amplitude A4{p —>■ p') is given by just one diagram with a shock wave crossing vertex 
(denoted by a cross): 

\x- P' 



M{p ^ p') = p I{p ,p — p') (p^ = p'^ = 0, p =p 



(4.5) 



As a more complicated example, let us compute one of the diagrams appearing in the compu- 
tation of the amplitude pi,p2 P3,Pi\ 



X 




We have two (p^ vertices and three shock wave crossings. The dotted lines stress the x~ ordering 
of the vertices. The fcj are the internal line momenta, whose — and _L components are fixed via 
momentum conservation, while the -|- components are determined by being on shell. There are 
two intermediate states: one before the shock wave {ki,k2,P2), and one after (^3,^3,^:4). The 
value of this diagram is thus: 



[ {2n)'6^'\J2 q« + Pi + P2 - P3 - P4) n 

^ a=l 



(2vr)^ 



X {-iXf k-^ I{k^ , qi) k^ I{k^ , q2) P2 HP2 > Qi) 
i 



X 



{Pt + Pt) - (kt + kt + pt) + IS (p+ + pt) - (p+ + k+ + kj 



n 



9{k- 



16 -^-^ fcj 

1=1 ' 



^In physical applications below we will be computing amplitudes with one incoming and several outgoing 
particles. This amplitude containing two incoming particles is considered for illustrative purposes only. 
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Finally, we have to discuss the relation between the transition amplitude M.{i ^ f) and the 
full relativistic scattering amplitude A^rei, i-e. the one obtained when we reinstate particle A as a 
quantum particle as opposed to replacing it with its classical field. We have: 

MU^ + 1 ^ A' + f) = -i2p\M{i ^ f) . (4.6) 

The incoming (outgoing) momenta of particle A are assigned as follows: 

Pa = (Pa'0'0)' Pa = (P^'O'Pi - P/) ■ 

In other words, particle A absorbs the total transverse momentum exchanged with the shock 
wave. As long as the momentum transfer is small compared to A' is almost on shell and the 
approximation is justified. 



The relative factor —i2p'\ in (4.6) is related to the normalization of the particle A state, 
which is lost when we replace it with a classical field. This factor is thus process- independent. 
For instance one can extract it from the external field approximation in QED [2^. The extra i 
can be traced back to the external field creation vertex, which carries a factor of i. 

This settles the question of relativistic normalization of the amplitudes computed via 't Hooft's 



method. We can now complete the comparison with the resummation method. Using Eqs. (4.5) 



(4.6) we have 

MUA + B^A + B') = -2tp\pl I{pl, q) , (4.7) 



which agrees with the eikonal amplitude from Eq. (2.2) including the normalization, modulo the 



difference between the S- and T-matrices already discussed in Section 2.2 



4.2 Gauge field 

In order to keep technical details to a minimum, we will not consider fermionic fields in the 
shock wave background. Instead, we will stick to a toy model in which charged matter (partonic 
constituents of the colliding hadrons) consists of massless scalars. This will be sufficient given our 
general goals. On the other hand, since the coupling constant of the 0^ lagrangian has dimension 
of mass, the cubic self-interaction is not a good model for the QCD radiation. We do have to 



introduce gauge fields. Thus we switch from (4.1) to a different microscopic lagrangian, describing 



the SU(3) Yang-Mills theory and a massless complex scalar in the fundamental representation, 
propagating in the shock wave background: 

^ = M\g^''{D^ct^yD,ct^ - ]^Ti[F^'' F^^]) . (4.8) 
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Most of the 0^ formalism is carried over with trivial changes. Instead of repeating the whole 
discussion, we will introduce the necessary modifications on a concrete example. 

Namely, let us consider the one-gluon emission: particle B, while scattering in the gravitational 
field of particle A, emits a gluon. The amplitude is given by the sum of the following two diagrams: 



P 




(I) 



P 



(11) 



p 



(4.9) 



The new objects are the gluon emission and the gluon shock wave crossing vertices. 

To simplify the computations, we will impose the Lorenz and light-cone gauge conditions: 



= 0, = , 



The treatment in a general gauge and demonstration of gauge invariance is given in Appendix |Aj 
Gluon emission in curved space is described by the cubic term in the lagrangian: 

(4.10) 



Here Qs is the strong coupling constant, and the SU(3) generators are normalized by Tr(T"T^) = 
1/2. In the light-cone gauge the singular component g^^ oc 5{x^) drops out (see Appendix [A| for 
a more detailed discussion). The gluon emission vertex is thus the same as in fiat space: 



P2,3 

Phi 



gsmPi+P2)-^ 



g,mp, + P2 - ^S^l).e (e+ = 0, l.e = 0) 



(4.11) 



where we used the Lorenz gauge to eliminate the £_ component. 

The gluon shock wave crossing vertex contains the same factor p^I{p^ , q) as in the scalar case. 
A new feature is that the £_ polarization component changes in the crossing according to: 



£2- 



£iq 

Pi 



, £2 = ei {e+ = 0) , 



(4.12) 
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This rule is easy to guess from consistency with the imposed gauge; see Appendix |X] for an exphcit 
derivation. Notice however that we don't have to keep track of this change in £_ if we use the 



simphfied gluon emission vertex in (4.11). 

We are now ready to evaluate the above two diagrams. Working for simplicity in the frame 
where p = 0, we get: 



(2p' + 1 



/(r,k)/(p'-,q-k)- 



fk-1 



l- 



[l-k)).e 



-[j9-/ + k]+-[/-k]- 



(4.13) 



Physical consequences of the derived expressions will be discussed below. 



As a final comment, we note that lagrangian (4.8) contains also a cubic gluon self-interaction 



vertex, which could be discussed analogously to (4.11), as well as two quartic vertices {(jxpAA and 
^4^74^). The quartic vertices do not contribute to the amplitudes in the coUinearly enhanced 
region, and we will not need their precise expressions. 



5 Initial state radiation in T-scattering 

The dominant QCD radiation effects in the usual perturbative hard scattering processes are the 
collinear initial and final state radiation. We now proceed to see how these effects manifest 
themselves in the T-scattering. We will focus on the initial state radiation and its effect on the 
parton distribution scale. Final state radiation, which happens after the partons cross the shock 
wave, is expected to be as usual. 

5.1 Observable 

To discuss radiative corrections to the PDFs, we need to choose a process and an observable which 
can be defined and computed beyond LO. The simplest such process is the T-scattering analogue 
of the DIS. In other words, we will consider an electron-proton T-scattering ep —>■ e-|-anything at 
a fixed momentum transfer. This is like in Fig. [TJwith an electron instead of a neutrino. 

The scattering is characterized by t = — and q~^, the energy transfer to the proton. These can 
be measured by observing the electron. As usual, we assume small angle scattering: \t\ <^ s. We 
will also assume that the relative electron energy loss is small, Pa- Under these conditions, 

and also since the electron does not QCD-radiates, we can represent it by a classical relativistic 
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point particle of fixed energy. This is our 'particle A\ Using the A' on shell condition {pa — q)"^ = 0, 
it is easy to show that the momentum transfer is mostly in the transverse plane, as expressed by 
the relation: 

q2 = q'{l + 2g+/p+) ~ _ 

Like in the DIS, we are interested in the differential cross section with respect to q and the 
Bjorken x: 

da 

X = -, < X < 1 . 



As is customary, we will first analyze the partonic cross section a between the electron and a 
quark (particle B). We will work in a toy model of scalar quarks. At LO (no gluon emission), the 



amplitude is (4.7) and the partonic cross section is given by 



5.2 One gluon emission in momentum space 

Armed with the formalism from Section |4], we can easily write down the gluon emission amplitudes. 



At the NLO we have diagrams with real gluon emission, as in Eq. (4.9 ), as well as virtual corrections 
to the external legs and the vertices in the elastic amplitude. As usual, the latter diagrams do not 
have to be computed explicitly, since they only correct the coefficient of 6{x — l)j^We thus focus 
on the real emission. 

The partonic cross section with one gluon emitted is given by a phase space integral (see 
Appendix [B| 

(^'5"NLO 1 f dH dz . 2/f - |2 



dz 

/ -^<5(x-qV(PB9+))|A^reir, (5.1) 

Jo ~ ^) 



d'^qdx IGn'^s J 2(27r)3 

g+ = iVr + (q - ir/Pn. • 
Here Airei is the relativistic scattering amplitude A + B A' + B' + g, related to the transition 



amplitude in the external field M.(B B' + g) via Eq. (4.6). The amplitude A4 is in turn the 
sum of the two diagrams (4.9), evaluated in Eq. ( |4.13 ). 



^Vertex corrections play a role in cancelling the IR divergence in the cross section corresponding to the emissions 
of soft gluons at large angles. This cancellation holds for any hard process. Specifically, it will also happen for the 
T-scattering because soft gluons do not feel the shock wave of particle A. Thus we do not discuss soft gluons in 
what follows, concentrating on the coUinear divergence. 
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We are using notation from (4.9) with p = pb, p' = Pb'- The q'^ is the total + momentum 
of the quark-gluon system after the colhsion. The z is the p~ momentum fraction carried off by 
quark B': 

p'~ = zp~ , = (1 — . 



Let us first analyze which region of the 1 plane contributes to the integral (5.1). For which 1 is 
there a z saturating the 5-function? The relevant function (see Fig. [?]) 

Xiz) ^ ci^ip-^q- 



(q - 1)7^ + 17(1 -z) 



has a maximum value 



max X(z) = Xiz^, 

0<z<l 



|q-l| + |l|)2 





q- 


-i| 




|q 


-1 


+ 


|i| 



Thus, the integrand of (5.1) is nonzero for I belonging to the ellipse: 

|q - 1| + |1| < |q| 



X . 



In other words, phase space limits the transverse momentum of the emitted gluon to be at most 
0(q). 




Figure 7: The function X[z) for a generic 1 (solid line) and for \\\ <^ |q| (dashed line). 



Let us now examine the amplitude, whose two parts are given in Eq. (4.13). By analogy 
with the usual DIS, we expect that part (I), corresponding to the gluon emission after the hard 
scattering, gives only a finite correction to the cross section, while part (II) contains a logarithmic 
IR divergence which has to be absorbed by redefining the PDFs. Let us see formally how this 
happens. 
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Notice that part (I) of the amphtude is non-singular in the 1 plane. In particular, the inter- 
mediate state denominator is completely fixed at q^(| — > 0. Omitting the x dependent 
factors, the amplitude is thus 0{gs /q.e/q^), and its square is 0{g'^\I\'^ /q^). After integrating over 
the ellipse in the 1 plane (area ~ vrq^), we get a finite contribution to the differential cross section 
of the relative order (9(as/vr). This is as expected. 

Interesting physics is associated with part (II), whose expression can be simplified as follows: 

= -igs^2p-ze.M, 
M. = J 1^ ~ /(^p-, q - k) /((I - z)p-, k) (5.2) 

+ / (^Mk-ip ^^"^"''^"''^^^^^"'^^''''^- ^^-^^ 

Here we separated the regular part / of the I{p~,q) from the 5-function piece describing the 
propagation without scattering: 

J(p-,q) = (27r)V2)(q) + J(p-,q), 

I{p~,q) = 27rblFM^\), 



where Fn{y) are the same functions as in Eq. (Q. We omitted a total 6^'^\q) piece from (^5.3|. 



The physical meaning of the decomposition (5.3) is as follows. In the first two terms, only one 



of the two splitting products of quark B participates in the gravitational interaction, the other 
one passing the shock wave without scattering. The last term instead describes their coherent 
gravitational scattering, as in Fig. [3j We call it the rescattering term, since it corresponds to the 
situation when the emitted QCD radiation changes its direction in the field of the shock wave. 

We now proceed to studying corrections to the cross section. Consider first the case |q| < b^^. 
In this case all the entering / functions are (9(27r6^). The rescattering term can be estimated by 
integrating up to |k| ~ beyond which point the / decrease faster than |k|~^, and the integral 
converges. We get 

|M.esc| ~ ^^e(2vr6^)^ ~ b^{2nbl) < ^^{2nbl) (|q| < b'^) . (5.4) 



We see that rescattering is subleading to the first two terms in (5.3). 
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Concentrating on the first two terms, the dominant contribution to the cross section comes 
from the singularities at 1 ^ and 1 ^ q. Squaring the amplitude and integrating we get: 

^ - i {|/(^r,q)|%^Q(^) + |/(xp-,q)r^Q^.(^)} |iog^ (|q| < 6-), 

Pq„.q(x) = Cp — , Pq^,(x) = Pq^q(1-x), Cf = AI?>. (5.5) 



Here we used that in the relevant regions of integration (see Fig. [7j) 

x = X{z)c^z (|l|<|q|), x = X(2)~l-z (|q-l|<|q|). 



Eq. (5.5) has the standard factorized form expected from an NLO QCD correction to a hard 



scattering [27]. The IR divergent logarithm multiplies the quark-electron and gluon-electron LO 



10 



cross sections, with the scalar quark DGLAP splitting functions Pq^q and Pq^g as coefficients 
As usual, we can absorb the IR divergence into the quark (first term) and gluon (second term) 
PDFs. If we fix the parton distribution scale at the upper cutoff, = q^, then the whole 
logarithmic correction is absorbed. 

It is of course not surprising that we managed to recover the standard factorization for |q| ^ 
b~^: rescattering was not important in this case, and without rescattering there is no difference 
between transplanckian and any other hard scattering. 

Let us proceed to the case |q| ^ b^^. The situation here is more complicated since the 
rescattering is no longer subleading. Consider for example the region |q| ^ |1| ^ 6^^. The 
rescattering integral is dominated by |k| < b~^, where J((l — z)p^,k.) is maximal, and not, say, 
by the region of |k| ~ |1|, where (^rjji is maximal. The reason is that /((I — z)p~,k.) decreases 
faster than Ikl^-*^ for |k| ^ 6^^. We get an estimate: 

|M.esc| ~ ^^(2vr&D \Iizp-,q)\ ~ ^|7(^p-,q)| (|q| » |1| » b^') , (5.6) 



which is comparable to the first term without rescattering in (5.3) 



5.3 Impact parameter picture 



In a situation when rescattering cannot be neglected, the separation into three terms in Eq. (5.3) 



becomes artificial, and we should treat the whole amplitude as given in (5.2). Substituting the 
definitions of /, we can transform this expression into a transverse plane integral: 

M, = — [ (fy'dPy ^-^' ~ y)v '-(q-')-y'+--V'^(yOe-W(i-^)V^(y) . (5.7) 

2vry |y'-yr 
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See [33] for the splitting functions of a colored scalar. 
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This simple equation provides a key insight into the physics of the process. Namely, we can view 
the factor \l/(y',y) = jy/Zyp as the coordinate-space wavefunction of the gluon-quark state into 
which quark B splits. Upon crossing the shock wave, this two-particle wavefunction is multiplied 
by the eikonal factors e*^^**^^'-* and e*'-^"^-'^**^^-', depending on the transverse plane position of 
each particle. Finally, to compute the S-matrix element, one takes the overlap with the outgoing 
state wavefunction e~''-^~'-''^'e~'''^l"1 



The completely general expression (5.7) can be further simplified if |1| <^ |q|. In this case, 
the typical y contributing to the integral are much larger than the typical y'. Approximating 
y' — y ~ — y, the amplitude takes a factorized form: 

M, ~/(^p-,g)f, (|l|<|q|), (5.8) 

where the gluon emission factor fj is given by 

^r = -^ I dV^e-^'-^+^(^-^)(^)" = _|l X C„[(l - zfl-hM , (5.9) 
ZTT J |y|^ F 

/■oo 







This factor has a nontrivial dependence on 1. For |1| ^ 6^ ^ integral is dominated by large 
|y| ^ so that the second term in the exponent can be dropped, giving 

f.^-Wi2 (|i| «&;'). 



Thus we recover the first term in Eq. (5.3), which in this limit dominates the amplitude. 

On the other hand, in the opposite limit the integral can be evaluated by stationary phase, 
with the result: 

pi X phase 

i,^'-^==h/\' (|1| (5.10) 
V + 1 

Remember that precisely in this case we expected a non-negligible contribution from rescattering. 



see Eq. (5.6). We now see that its effect is indeed important: rescattering leads to an 0{1) 
reduction of the gluon emission amplitude! Physically, this can be explained as follows. To 
scatter with large 1, the emitted gluon must cross the shock wave in the region of small impact 



parameters. In this region, the eikonal factor in (5.9) distorts the gluon wavefunction, which leads 



to suppression of the amplitude via destructive interference. It is however remarkable that, up to 
an 1-dependent phase, the suppressed amplitude still goes as Ij/P. 



^^See |34j where similar considerations axe made for a gluon splitting into a qq pair in the field of an incoming 
nucleus. 
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In terms of the function C„, the above asymptotics can be stated as follows: 
Cniu)-^! (m^O), \Cn{u)\ ~ 1/Vn + l (m>1). 
As can be seen from Fig. Isl the large u behavior sets in already for u > 2. 
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Figure 8: Numerical plots of Cn{u) for n = 2,3,4,5: the absolute value (solid blue 
curve), the real and imaginary parts (dashed and dot-dashed curves), the u = oo asymp- 
totic value l/y/n + 1 (dotted line). 

We are now ready to derive the logarithmic correction to the cross section for |q| ^ b^^. We 



square Eq. (|5.8|) and integrate in |1| <^ |q|, taking into account the suppression for |1| ^ b^ . 

'I'll dl 



Schematically, we get: 



dH\i\' 



dl 1 

+ 



log 



/^F(|q|) 



(5.11) 



'|l|<|q| ^WR ^ I /XiR 

In other words, due to the l/(?7, -|- 1) suppression of the second term, the arguments of the two 



logs combine into a geometric mean which coincides exactly with Eq. (1.1)! 



The final result is as follows: for |q| ^ 6^ ^, the NLO correction to the cross section is given 



by the same Eq. (5.5) as for |q| ^ &c ^ with the following simple replacement: 

log^^ (|q|»r' 
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5.4 Discussion 



Let us discuss what we have achieved. First, we have shown that the emission amphtude of near- 



colhnear gluons has a factorized form, Eq. (5.8). Physically, it means that we can first consider 
the hard scattering process, and worry about adding additional gluons later. If a parton splits in 
two before crossing the shock wave, one and only one of the splitting products absorbs most of 
the momentum transfer in a typical event. This is actually an important check of validity of the 
partonic picture. 

Second, we found explicitly the gluon emission factor fj. The probability distribution of emitted 
gluons is given by |f(p. We found that for large relative transverse momenta |1| ^ (but still 
|1| ^ |q|) this distribution is suppressed by a factor l/{n + 1) relative to the standard QCD 
distribution This is a new effect, which could be used to measure the number of extra 

dimensions. 

Finally, as a consequence of this suppression, the logarithmic NLO correction to the partonic 
cross section involves, for |q| ^ b~^, a scale which interpolates between the usual |q| and in 



agreement with Eq. (1.1 ). In fact, as we show in Appendix [C| such logarithms occur in every order 
of perturbation theory. Thus, as usual, they can be exponentiated and removed by shifting the 
parton distribution scale to /Uir(|q|). This, then, provides a formal justification for the proposal of 
Emparan, Masip, and Rattazzi [17], that this scale is the one minimizing higher-order corrections. 

Our derivation of the scale (1.1) has an added advantage that we now know the gluon distri- 
bution. This distribution could not be easily guessed: for example the standard 1/q^ with a sharp 

cutoff would give rise to the same log. At the same time, the very fact that we found agreement 
with Ref. [IT] may seem like a mistery. Remember that they fixed this scale to be equal to b^^ 
(see Section [3]). Where, then, does the typical impact parameter b^ hide in our computation? 

In fact, if we are not interested in the gluon distribution, we can reformulate our derivation 
so that it will conform with the original intuition of [17J. The idea is to compute the LHS of Eq. 



(5.11 ) not from the asymptotics of Cn but directly from the definition of fj in the impact parameter 
representation. Since the L2 norm of the gluon wavefunction is the same in the momentum and 
in the position space, we have 



I'^l dH 



(27r) 



dy , y 
oc log 



max 



The only subtlety is that in the LHS we are not integrating over the whole 1 plane, and thus the 
limits of the y intergration have to be adjusted accordingly. A moment's thought shows that the 
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correct limits should be put at the typical y values contributing to fj at |1| ~ //jr and |1| ~ |q|: 

Umax ~ A'lR, ) Umax ~ ^* ■ 

So, we recover the same logarithm as above, and this time b^, enters explicitly. 

We have worked throughout in the toy model of scalar quarks. However, it should be easy to 
adapt our considerations to the realistic case of fermionic matter. One would have to compute 
the shock wave crossing vertex for the fermion field. This will require solving the Dirac equation 
in the shock wave background. We expect that our results about factorization and suppression of 
radiation at large angles will remain true in the fermionic case as well. 

6 Simultaneous radiation 

So far we were making the technical simplifying assumption that particle A does not QCD-radiate. 
This is of course not true in pp collisions at the LHC, when both coUiding partons are colored. 
We will now discuss how one could relax or remove this restriction. 
Consider the following two key properties of QCD radiation: 

1. Near-coUinear emission dominates. 

2. Its amplitude takes a factorized form. 

These properties are true for the standard hard perturbative processes. For the T-scattering 
with non-radiating A, we have also shown them to be true, provided that a gluon emission factor 
is adjusted accordingly. We conjecture that these properties continue to hold when both A and B 
are allowed to radiate. In practice, this implies that the dominant part of the emitted radiation 
can be described using the independent emission approximation: first compute the radiation off 
B (taking A classical), then off A (taking B classical). 

Intuitively, this can be justified as follows. The fact that the near-collinear emission dominates 
is due to the ~ 1/|1| singularity in the gluon emission factor, combined with a phase space cutoff 
\M ~ IqI- We have seen that shock wave crossing tends only to suppress radiation at large angles 
by distorting the gluon wavefunction, and we do not expect this tendency to reverse when A is 
allowed to radiate. Thus collinear radiation should dominate also for the T-scattering. Once we 
know that collinear emission dominates, and thus the emitted radiation does not change the hard 
transverse momentum flow of the process, it seems reasonable that the amplitude should factorize 



28 



into the product of the 2 — 2 scattering with the hardest momentum exchange times the gluon 
emission factors!^ 

The independent emission approximation is probably adequate for most practical purposes. 
Nevertheless, to try to go beyond it is an interesting theoretical challenge. We will now describe 
pictorially a generalization of our formalism which, we believe, can describe simultaneous emissions 
from A and B without any extra approximation (except, of course, large CM energy and small 
scattering angle). 



The starting point is the emission amplitude in the impact parameter representation (5.7) 
(see Appendix [C] for its generalizations to two and more gluons). The physical meaning of this 
equation in terms of the two-particle wavefunction and individual eikonal factors was explained 
above. Suppose now that both particles split. The amplitude can be constructed according to the 
following three steps (see Fig. [9]). 

1. We evolve the partons from infinity to the transverse plane x~ = x'^ = where the collision 
is assumed to happen. We introduce many-particle wavefunctions for the splitting products 
of both A and B. The total wavefunction in the transverse collision plane is the product of 
the two: 

^toi = ^^({xJ)^B({y6}) . 

Here (y^) are the transverse coordinates of left- and right-movers. The can be 

computed via the flat-space light cone perturbation theory. For the one-gluon emission they 



are the same as in (5.7). 



2. When the splitting products cross the transverse collision plane, the wavefunction "^tot is 
multiplied with eikonal factors, one for each pair of opposite-movers. It is not difficult to 
guess that these factors are equal to 

expiZaZbXi^a - Yb) , 



where x is the 2 — > 2 eikonal phase from (2.3), and Za = Pa Ip\i = Pb IPb 



longitudinal momentum fractions carried by partons a and h. 



^^We are ignoring here soft gluons which may be exchanged between particles A and B. In case of standard 
haxd QCD processes involving two initial hadrons, like in Drell-Yan, where proofs of factorization to all orders in 
perturbation theory exist [35 , it is known that such soft gluon exchanges cancel in the total cross section. We 
expect this cancellation to carry over to our case, because soft gluons do not feel the gravitational field of the 
energetic particles (the eikonal factor being proportional to the gluon energy). 
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3. Finally, to find S-matrix elements, we compute the overlap with the outgoing state wave- 
functions. These are simple plane waves if the partons do not undergo any splittings after 



the shock wave crossings. If such splittings are present, like in part (I) of (4.9), these are 
more complicated functions of external momenta. However, since at this stage left- and 
right-movers no longer interact, the flat-space light cone perturbation theory can be used to 
find them. 




Figure 9: This diagram represents one of the terms in the amplitude A + B ^ A' -\- 
g + B' + g, in which both the right-moving A and left-moving B split before colliding. 

We think that this generalized formalism, apart from providing an attractive mental picture, 
could find interesting future applications, especially in the problem of gravitational radiation 
emission. 



7 Conclusions 

The main point of this paper is that including radiation in small-angle transplanckian scattering 
is, after all, a tractable problem. The gravitational interaction producing eikonal phase factors 
happens instantaneously, while the processes of particle splitting are spread in time. Based on this 
observation, we developed a formalism which allows explicit computations of scattering amplitudes 
in presence of hard, quantum, radiation. In the impact parameter representation, the radiating 
particle splits into a multi-particle virtual state, whose wavefunction, computed via light-cone 
perturbation theory, is then multiplied by individual eikonal factors. 

We demonstrated the usefulness of our formalism on the concrete problem of initial state 
QCD radiation in transplanckian scattering. We will not repeat here the detailed discussion of 



the obtained results given in the Introduction and in Section 5.4 We believe however that this 



example by no means exhausts the list of possible applications. We are particularly hopeful about 
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a possibility to shed new light on the problem of gravitational radiation emission, which is always 
on our mind. 
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A Gluon field in the shock wave background 

Here we will discuss quantization of the free gluon field in the shock wave background, as well 
as the gluon emission vertices, the gauge invariance, and the absence of gluon emission terms 
localized on the shock wave. 

We start from the Maxwell Lagrangian in curved spacetime, which integrating by parts can be 
rewritten as 

Notice that while the D-dimensional Ricci tensor of the AS metric is zero for x_l 7^ as a 
consequence of Einstein's equation, the 4-dimensional Ricci tensor of the metric g^y induced on 
the SM brane is nonzero. Namely, it has a nonzero component 



(compare with (2.6) where the Laplacian is with respect to all D — 2 transverse directions). 



Let us fix the curved space Lorenz gauge D^A^ = 0. In this gauge the EOM take the form: 

D2a,-(%/A^ = 0. (A.l) 

We first discuss the gluon propagator across the shock wave. To find it we need to solve the 
EOM with the initial conditions = e^e*^^ for < 0. 

Consider first the light-cone gauge case = 0. It is easy to see that the shock wave metric 
has r^l^j, = 0, which implies ^4+ = 0. The other gauge field components are found to satisfy the 
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equations: 



4 



Just as in the scalar field case, it is enough to find a matching condition, i.e. to relate solution 
of these equations for x~ = — e and x~ = +e. One could use the method of Section 2^ based on 
transforming to the x coordinates. Here we want to demonstrate a different, equivalent, approach. 
Namely, we regularize the above equations by smearing the 5-functions. One can show that the 
terms put in the RHS of the equations can be dropped in this analysis, since their effect goes to 
zero when the regulator is removed. All the other terms are however important. 

From the first equation, we find how the transverse components varies across the shock wave: 



e e 



hip.x 



exp(z- 



5(x-)$(x)). 



(A.2) 



Subsituting this solution into the equation for A_ we find: 

1 



A_ = [e_--e{x-)eAH^)]e 
Thus on the other side of the shock wave: 

A(x~ = +e 

A_{x' = +e] 



-i^^x'^ +ip.x 



(5(x-)<l>(x)). 



(A.3) 



£ ^-i'^(x+-'^(x))+ip.x 



— $(x))+ip.x 



These are the desired matching conditions. Taking the Fourier transform we recover the rule 



(4.12). 



A general solution of the EOM (A.l) is a linear combination of the just found solution in the 
gauge A^ = with a pure gauge solution 



An = d^ip . 



(A.4) 



The gauge parameter ip is given by Eq. (4.2) as a general solution to the Klein-Gordon equation. 



Let us now discuss the gluon emission term (4.10) in a general gauge. There are several 



quantities in (4.10) with 5- function singularities on the shock wave. Thus one may wonder if there 
is a contact emission term localized on the shock wave. In fact such a contribution is absent, so 
that one can always compute the gluon emission as a sum of two separate integrals for > 
and a;~ < 0. To see this, one can argue as follows. 
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First of all, as already mentioned in Section 4.2, localized terms are absent in the light-cone 
gauge y4+ = 0. In this gauge does not contain 5-function singularities as one can see from 
the explicit solution (A.2),(A.3). The 5-function does appear in g~^^ and in (p* d . 



but all 



contractions involving these terms necessarily contain and vanish. 



Second, consider the pure gauge case (A.4). In this case there are several 5-function terms in 



(4.10 ). However, one can show that they cancel among themselves. The reason for this cancellation 



is as follows. Since the integrand in (4.10) is a Lorentz invariant, we can compute it in the x 



coordinates (2.8). In these coordinates both and %jj are continuous, and the integrand has at 



most 6'-function singularity on the shock wave. 

To demonstrate the absence of localized terms by a concrete example, let us show that the 
one-gluon emission amplitude is gauge invariant. We thus have to show that the amplitude to 
emit a longitudinally polarized gluon is zero, without inclusion of any extra terms localized on the 
shock wave. The one-gluon emission amplitude is given by the coordinate-space integral: 



M = z d'x^g'^''{[r/ix)Y^,4>f{x)]A:-\l,s-x). 



Here A°^^(l, e; x) is the outgoing gluon wavefunction. In the considered longitudinal case ^ ^ — 



lOUt , 



X) 



we have: 

A-*(/,£;x)(x V« 

The integral splits into two parts: x^ > 0, x^ < 0. Each of these can be integrated by parts and, 
using the current conservation, reduces to a boundary term localized on the shock wave. These 
boundary terms are not quite identical: 



M(x->o) 

-^(x-<0) 



dx+d^x <i 



^t[]^p'-x+-p'.x) Q g-i(|p-[a;+-$(x)]-p.x) 1 i(iZ-x+-l.x) 



j dx+d'^ |g<|p'-[x++*(x)]-p'.x)^ 



g-*(|p-x+-p.x) 1 gi(l/-[x++$(x)]-l.x) 



However, after integrating in x~^ and taking into account the resulting p -conserving (5-function, 
they are seen to cancel. 



B Light-cone phase space 

The partonic cross section with n gluons emitted A + B ^ A' -\- B' + gi + . . . + Qn is given by the 
phase space integral 
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We assume that A and B collide head on along the z direction. We use the light-cone phase space 
adapted to the direction of motion of each particle: 

-P/) = 25(p+ -j9+)5(p,r -pj)5(2)(p^ _ p^) . 
The momentum conserving (5-function is saturated by integrating in dp\id'^'pA'dp']^,., which gives 

1 An), 2 Pb i2 "TT dlTdHi 



The differential cross section in the momentum transfer q and the Bjorken x is thus given by: 

da _ 1 /" I A^(")|2 TT dl^dHi f q2 



d\dx IQir^P J ' p^, J-1 (27r)32/- V p-q^ 

where g+ = p\:/p~b' + Pb' = Q - li, P^' =Pb~Y1 K under the integral sign. 

C Multi-gluon emission 

Let us look at the leading logarithmic corrections to the cross section which come from the radiation 
of many gluons. We want to verify that they have the correct form in order to be reabsorbed into 
the PDFs normalized at the scale fipiQ)- 

To begin with, consider the emission of two gluons. In the LLA, we are looking for (a^ log 
corrections to the cross section. One can show that, because of the form of the denominators in 
our Feynman rules, a gluon which does not cross the shock wave will not give rise to a large 
logarithm. Thus the relevant diagrams are those with both gluons emitted at x~ < 0, shown in 



Fig. 10 Let us consider the first of these diagrams, and define: 



P'~ = ZiZ2P~ 

= (1 - Z2)p~ 
= Z2{1 - Zi)p~ . 

Computing this amplitude using the Feynman rules and going into the impact parameter repre- 
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q + JI; k; - JI; 




q + Eik; -E; li 




Figure 10: Relevant diagrams for the emission of 2 gluons. 
sentation, we find A^^'^'i) _ 



Bp'-(2^?,)2Mj'^)ei,e2,-, with: 



M 



{hh) 

ij 



:T^ki+k2) 



(27r)2 {2ny {k^y (^k^^ + '-^f^k^' + 2k, ■ k^) 

X I{ziZ2P~,q + ki + k2 - li - I2) /(2;2(1 - zi)p',\i - ki) /((I - Z2)p' - k2) 
1 /" .0 .0 .0 fx-vi)i (x - y2)j + i?(x - yi)j 



(27r) 



r2 ,2 ,2 x-yii 

d xrf yid y2-^ 

x-yir 



X e 



(x-y2) + 5(x-yi)|2 + A|x-yi|2 

-i(q-li-l2).x+i2;i22^"3?(x) -ili.yi+i2i(l-Z2)^'I'(yi)p-il2-y2+«(l-2:2)^"3?(y2) 



where A = zi{l — Zi)/{1 — Z2), B = 1 — zi. The integrand has the expected form (three-particle 
wavefunction in the transverse plane) x (individual eikonal factors) x (outgoing states). 

Let us estimate this amplitude for |q| 3> |li|, |l2|- In this limit we will find the double logarithm 
associated with Pq^q{zi)Pq^q{z2), while in the other regions there are those associated with 
splittings into gluons. Now, the integral in x is dominated by values much smaller than those 
dominating the integrals in yi and y2. Neglecting x with respect to yi and y2, we obtain: 



I{ziZ2P ,q 
1 



1 



(27r)2 



y2j + Byij 



|yi 



(y2 + 5yi)2 + Ayi2 



-d2y2 i(l-^2)(^^ 



A double logarithm can be obtained only if the amplitude behaves like |M 



this means that we must look for the situation in which the integrand behaves like 



1 1 
|iil II2I 

|yiP Iy2| 



(C.l) 



. In turn, 



2 ) 



which 



35 



is true only if the yi integral is dominated by much smaller values than those dominating the y2 
integral. This will happen if and only if |li| ^ |l2|. We conclude that the gluon emitted closer to 
the shock wave must have a larger transverse momentum. We will come back to this later, when 
we generalize to an arbitrary number of gluons. 

Thus for I 111 ^ II2I, we can neglect yi with respect to y2 in the integrand, and the last integral 
factorizes giving: 



M 



(hh) 



I{ziZ2P ,q) fi(/i ,li) ij{l2 ,h] 



|ll|»|l2|) 



where fj is the gluon emission factor introduced in Eq. (5.9). After squaring and integrating in li 
and I2, this gives a double logarithm of ^f{(i)- 

For completeness let us verify explicitly that the contribution from the region |li| ^ II2I is 
subdominant. One can show that the amplitude depends on li and I2 as follows: 



m£^'^) oc <; 



l2jl2i 



|l2|4 



FT I ^ 
II2I 



1+— !— / 1^ 

log / I'^l^ 



I 111 be' 



llll* VI12I 



if |li| < II2I < 
if |li| < b-^ < II2I 

if b-^ < llil < II2I 



It is immediate to see that from each of the three regions we only get single logarithms. 
Thus, the region II2I < |li| < |q| dominates, and in the LLA we get: 




d%d%\M 



ij I 



1 



I{ziZ2P ,q)| < vrlog 



where the ^ factor comes from the — II2I). 



The second diagram in Fig. [TO] is the same with li ^ l2- Since the leading contributions come 
from integration over different regions of transverse momenta, in the LLA there is no interference 
between the two diagrams. Putting all together and including also the double logarithms arising 
from the regions |q| ^ |q — li|, II2I and |q| ^ |li|, |q — I2I, which can be computed in a similar 
way, we finally obtain: 

dd' NLO 1 



d'^q^dzidz2 47r2 



\l{ZiZ2P ,q)\'^PQ^Q{zi)PQ^Q{z2) + \l{{l-Z2)p ,q)\'^PQ^g{l - Z2)Pq-.q{zi] 



+ \I{Z2{1 - z,)p~,qfPQ^,il - Z^)PQ^Q{Z2)] ^ 



;^log — ^ 



1 2 



which shows that the optimal factorization scale is fipiQ)- 
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Figure 11: Typical relevant diagram for the emission of m gluons. 
Generalization to an arbitrary number m of emitted gluons is straightforward. Consider the 



typical relevant diagram shown in Fig. 11, and let us look at the transverse momenta involved in 

the quark-gluon vertices. In the first one from below only km appears, thus the contribution to 

the amplitude will be 1^ , like in the case of one gluon emission. In the next vertex both km 

and km-i are involved, but it is obvious that its contribution should reduce to ^TJ^ ^^"12"^ when 

1*^111— 1 1 

I km I ^ |km-i|- This is true for all the m vertices: if the transverse momentum k of a gluon is 
much larger than those of the gluons which were emitted previously, then its contribution to the 
amplitude will be Thus in general we can write, for the diagram in Fig. [TTj 

m 
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where: 



fi (2^)7 Vii + E.>, (^(kj ■ kr) + E.,.>, (^(ks ■ kr 



X 



m 



j=i Vj=i 

Fourier-transforming this expression and considering the hmit |q| ^ |lj| (j = l,...,m), we would 



obtain an expression analogous to (C.l) with m integrals instead of two. From the structure of 



the amplitude in k space, it is clear that it will be: 



where the multi-particle wavefunction 



^(yi,-,ym) ^ n^N^ if |yi| < |y2| < ... < |yr 

.7 = 1 ^-^j' 



This means that the only term with m^^ power of a large logarithm comes from the region |li| ^ 
II2I ^ ... ^ |lm|. When computing the total cross section with m identical gluons in the final 
state, we can always reorder the gluons so that |li| > II2I > ... > |lm|. Then in the LLA only the 
shown diagram contributes, with a ^ factor from 6'(|lciose| > ... > |lfar|). Putting all together, we 
get 

7-^ -1 ™ ™ -ir 9/\nm 

"CTnLO 



2tt ^'^ 



^ j=i i=i 

where: 

m 

i=j 

All these large logarithms are absorbed into the PDF normalized at the scale /ii?(|q|). 
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